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Ionic-liquid gates have a high carrier density due to their atomically thin electric double layer (EDL) and 
extremely large geometrical capacitance C g . However, a high carrier density in graphene has not been 
achieved even with ionic-liquid gates because the EDL capacitance C ED l between the ionic liquid and 
graphene involves the series connection of C g and the quantum capacitance C q , which is proportional to the 
density of states. We investigated the variables that determine C ED l at the molecular level by varying the 
number of graphene layers n and thereby optimising C q . The C EDL value is governed by C q at n < 4, and by 
C g at n > 4. This transition with n indicates a composite nature for C ED l- Our finding clarifies a universal 
principle that determines capacitance on a microscopic scale, and provides nanotechnological perspectives 
on charge accumulation and energy storage using an ultimately thin capacitor. 

Single-layer graphene (SLG) has two kinds of singularities. One is a topological singularity at K points in 
non-doped SLG, called the Dirac point. The carrier transport at the Dirac point, with a vanishing density of 
states D(E), has been intensively studied 15 . The other is the van Hove singularity at M points in highly 
doped SLG, which has recently attracted significant attention due to predictions of several novel phenomena such 
as superconductivity 6,7 , ferromagnetism 8 , and charge/spin density waves 9,10 caused by divergent D(E) and a band 
structure with high symmetry. We can reach this state by doping SLG with 0.25 electrons or holes per carbon atom 
(~ 1 X 10 15 /cm 2 ). This high carrier density is expected to be attainable with an ionic -liquid gate because of its thin 
electronic double layer (EDL) 11 " 14 , i.e. large geometrical capacitance determined by C g = Si L s 0 /d. Here, g IL , e 0 andd 
are the relative permittivity of the ionic liquid, vacuum permittivity, and thickness of EDL, respectively. However, 
studies have demonstrated that carrier doping into SLG with ionic liquids is much lower than expected 15,16 . This is 
because the potential difference between an ionic liquid and SLG is made smaller than the applied gate voltage V g 
by the shift of the Fermi energy £ F when doping SLG with an additional charge of Q. This phenomenon can be 
interpreted in terms of the quantum capacitance C q , defined as C q = eQAEp = 6 2 D(£f) 15,17,ls , which is connected to 
Cg in series. Thus, the total capacitance C EDL between the ionic liquid and SLG is understood as follows: 

1/C E DL = 1/C g +1/C q 



(1) 



The effect of C q is prominent when C q <C C g , which is characteristic of the combination of SLG with small C q and 
an ionic liquid with a large C g . Thus, even if an ionic liquid has a large C g , small C q limits C E dl to make high charge 
density unattainable. Here, increasing the layer number n might solve this problem, because D{E) increases with 
n, and accordingly C q increases. To find an optimal n for carrier doping with an ionic liquid is thus quite 
important for the discovery of novel ordered states in graphene. In addition, the C EDL will increase with n, but 
the specific surface area (the surface area per unit of mass) of few-layer graphene (FLG) decreases. In a practical 
viewpoint to realise the superior characteristics in graphene-based supercapacitor 19,20 , therefore, it is of signifi- 
cance to elucidate the optimal n producing both the large C E dl and specific surface area. In this letter, we examine 
the capacitance of SLG and FLG while systematically changing n. 

Results 

Our results are summarized in Figs, la and lb. The capacitance C E dl between the ionic liquid and FLG was 
measured as a function of n by two methods, a transport measurement using graphene field effect transistors 
(FETs) (Fig. la), and a capacitance measurement using AC (Fig. lb). When the gate voltage Vg induces the charge 
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Figure 1 | Layer-number dependence of EDL capacitance between graphene and ionic liquid (a), C EDL evaluated from conductivity measurements (Fig. 2 
and Supplementary Fig. S2 online) for SLG and FLG devices. The C EDL values measured in hole-carrier and electron-carrier regimes using two-terminal 
devices are indicated by blue solid circles and red open circles, respectively. The C EDL values evaluated from Hall coefficients are indicated by green 
triangles, (b), C EDL directly measured with AC (Fig. 3). The minimum C EDL values from V tg = — 1 to 1 V are plotted. In (a) and (b), the error bars are the 
standard deviation of the data measured for each n device. The black solid lines in (a) and (b) represent calculated C EDL based on a simple theoretical 
model (Fig. 4). Dependencies of experimental data on n are consistent with the model. 
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Q(V g ), C EDL DC = Q/V g and C EDL AC = dQ/dV g are obtained by the 
transport and capacitance measurements, respectively. These two 
measurements are complementary to each other and are especially 
important in the estimation of C EDL , which depends significantly on 
measurement frequency due to the slow transport of the ions in an 
ionic liquid. The calculated values of Cedl are also shown in Figs, la 
and lb, and are in good agreement with the experimental data. At 
n = 1, a small Cedl is observed, as in previous studies 15,16 . With 
increasing n, Cedl increases and saturates at n > 4. In the following, 
we describe the experimental and theoretical details. 

First, we evaluated Cedl from transport measurements. Two-ter- 
minal graphene FETs were prepared on Si0 2 300 nm/highly doped 
Si substrate. The structure of the device and the optical microscope 
image are shown in Fig. 2a. In order to obtain a reliable n depend- 
ence, we used only uniform graphene layers without any overlaps of 
different layers, wrinkles or folds. If the device has overlaps of dif- 
ferent layers, the exposure of interstices between the neighbouring 
graphene layers to the ionic liquid might lead to inaccurate n and 
errors in estimation of surface area. The uniformity of graphene layer 
was clearly confirmed by the optical contrast of micrographs 16,21-24 . 
Raman spectroscopy mapping and atomic force microscope (AFM) 
image also evidenced the uniformity of the graphene layer. The 
number of devices used in this measurement is shown in Table 1; 
good reproducibility was obtained. The sheet conductivity a was 
measured as a function of back gate voltage V bg applied to the Si 
substrate. Next, a droplet of ionic liquid was placed on the graphene 
surface, and a was measured again as a function of the top gate 
voltage Vtg applied to the ionic liquid with a Pt wire (100 urn in 
diameter). In Figs. 2b-2d, ffs of FLG FETs with n = 1, 4, and 9 are 
plotted as functions of V bg and V tg . Note that the horizontal scales are 
different for V bg and V tg . It was found that a was modulated with a 
lower gate voltage for V tg than V bg . The conductivity varied with gate 
voltage in proportion to the capacitance as follows, \da/dV t ^\ = 
CEDLrtg or | rfcf / ciVbg | = Cbg^b g > where /.( is electric field mobility 
and the subscript tg (bg) means top (back) gate. Thus, Cedl is given 
by 

da I d V tg /< bR 



C, 



■C bg 



(2) 



Here, <io"/ciV t g( b g)is defined as the steepest slope of the experimental 
data around the charge neutrality point. Assuming /i b g//'tg = 1 and 
substituting the experimental value of C bg = 11.5 nF/cm 2 , Cedl is 
calculated for each device and plotted in Fig. la. A clear n depend- 
ence of Cedl is found in Fig. la. Cedl is a minimum at n = 1. With 
increasing n, C E dl shows a maximum at n = 4 and decreases 



smoothly as n approaches 10. Maximum carrier density was esti- 
mated to be 0.003 /C atom (4.7 X 10 13 /cm 2 ) at n = 4 and V tg = 2 V. 

In order to prove the validity of Cedl(«) obtained by the transport 
technique, we directly measured the capacitance using AC. The 
device structure and the optical microscope image are shown in 
Fig. 3a. Uniform graphene layers are carefully chosen in this mea- 
surement. The graphene layer was covered with photoresist except 
for a region about 15 urn square, on which a droplet of the ionic 
liquid was placed. AC voltage was applied across the interface 
between the ionic liquid and FLG, and the capacitance Cedl was 
evaluated by the out-of-phase components of the detected current. 
Cedl was measured with an excitation amplitude of 10 mV at a 
frequency of 500 Hz and with a superimposed DC bias voltage V, g 
ranging from — 1 to + 1 V. The equivalent circuit of our system is 
shown in Fig. 3a 16 . Here, it should be noticed that the contact area 
between the ionic liquid and Pt electrode is much larger than the area 
between the ionic liquid and FLG. Hence, the impedance of Cp t is 
negligible compared to that of Cedl- This means that V tg is effectively 
applied at the interface between the ionic liquid and FLG. Note that 
V, g is uniform along this interface, because the impedance of C E dl (~ 
10 s Q.) is much larger than the resistance of FLG (10 3 - 10 4 Q.). The 
measured impedance is typically Z ~ 10 6 - 10 s i Q at a frequency of 
500 Hz. Therefore, the capacitive component dominates the total 
impedance, showing that Cedl can be directly determined by the 
measurement with LCR meter as described in the method section. 

In Fig. 3b, Cedl for fx = 6 is plotted as a function of V tg . We found 
that C E dl s for all the devices significantly depended on V tg , which 
corresponds to the ambipolar behaviour in D(E). Here it is important 
to notice that V, g in our experiment may be different from the actual 
voltage applied to the FLG since the potential drop caused by the 
interface of the Pt wire is not considered. The accurate dc bias in the 
electrochemical impedance analysis may be defined as the potential 
difference from the reference electrode using a three-terminal 
setup 15 . However, in our devices used in this study, the potential drop 
at the Pt interface is expected to be so small that applied voltage from 
the ground is approximately equal to the actual voltage, because the 
Pt electrode is much wider than the channel region of FLG, as 
described above. To compare with the transport measurement, the 
Cedl for each device, obtained at the charge neutrality point, is 
plotted as a function of n in Fig. lb. Here the charge neutrality point 
is defined as the V tg where Cedl exhibited the lowest value, and is 
located around V tg = 0 V as shown in Fig. 3b. Since V tg ~ 0 V, the 
value of minimum capacitance is obtained without any ambiguities 
of the bias voltage. Consequently, the discussion on n dependence of 
capacitance shown in Fig. lb is appropriate. As shown in Figs, la and 
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Figure 2 | C EDL obtained by 2-terminal transport measurement (a), Schematic side view of two-terminal graphene FET and optical microscope image of a 
device. Scale bar is 10 (im. Materials and terminal configurations are indicated. After <x( Vbg) was measured using the circuit shown by orange lines, an 
ionic liquid was applied and a{ V tg ) was measured using the circuit shown by blue lines. Actual size of Pt wire (100 urn in diameter) is much larger than the 
channel length (10 um), but it is depicted as a thin wire in the figure for clarity, (b), (c), (d), Comparison of a( V"b g ) (orange line) and a( V tg ) (blue line) for 
SLG (b), 4-layer (c), and 9-layer graphene (d). Note that the horizontal scales are different for V"bg ana V tg . To determine C EDL from equation (2), 
we estimated da/d V tg( b g ) as the steepest slope of the a( V tg (b g )) curve. The gate voltages were swept from negative to positive values. 



lb, a C EDL - n plot is substantially consistent with that from the 
transport measurement, which supports the validity of our estima- 
tion of C EDL from the transport measurement. 

Discussion 

Here, we discuss the validity of the Cedl obtained by two-terminal 
transport measurements. The analysis using equation (2) might 
oversimplify the estimation of Cedl because the contribution of 
fibgUkg is bypassed. First, we consider the validity of the assumption 
of /4,g/^tg = 1- This assumption is required because C EDL and ji tg 
cannot be determined independently without the result of Hall effect. 
In order to investigate the validity of this assumption, we carried out 
additional transport measurement using the multi-terminal Hall bar 
devices with n = 1, 2, 7, and 12. (see Supplementary Fig. S2 online). 
In this method, we can independently determine Cedl an d ^tg 
without any assumption. By measuring the Hall coefficient, ft H > 
and a in four-terminal geometry under the gate voltage V g (V bg 
and V t g)> we can evaluate j.i ((j bg and /( tg ) and C (Cb g and C EDL ) 
independently by using the following equations 25-26 . 

lt=\trRa\ (3) 

CV g =-l/« H (4) 

These equations are satisfied in the doped region where only one 
kind of carrier presents. i? H was measured as a function of V tg under 
perpendicular magnetic field of 6 kOe. The values of Cedl were 
estimated from the slope of IARh - V, g curve in the range where 



equation (4) holds. The C ED l - " plot is shown in Fig. S2b and 
C E dl s are added in the graph in two-terminal measurement 
(Fig. la). This shows a good agreement between the C E dl values 
obtained from two-terminal transport with assumption of /i bg //(tg 
= 1 and Hall effect without any assumption. Therefore, this result 
proves the validity of the assumption of ^bg/A'tg ~ 1 in the analysis in 
two-terminal geometry. 

Furthermore, we deduced ji^ and ,u tg separately by substituting 
both Rh and a in four-terminal geometry into equation (3), and 
evaluated i(b g /|U tg . The result is shown in Fig. S2c. The /x b g//'tg scarcely 
depends on n and is regarded as a constant around 0.5. The differ- 
ence, 1 and 0.5, in /( bg / i( tg between two analyses described above may 
be closely associated with the contact resistance. The contact resist- 
ance reduces the conductance measured in two-terminal geometry, 
and this effect becomes more significant in highly doped graphene. 
Accordingly, /.i t g in two-terminal measurement could be underesti- 
mated, because ii tg is estimated in more carrier-accumulated region 
than ii bg , while f.i tf , and ^ b g evaluated from _R H and a measured in 
four-terminal geometry are not affected by the contact resistance. 



Table 1 Num 


ber of FLG devices prepared 


in this study 






layer number 


1 2 3 4 5 6 


7 8 9 


10 


total 


Transport 
measurement 
Capacitance 
measurement 


4 4 4 3 1 5 
12 14 2 3 


1 2 3 
1 2 0 


3 
2 


30 
18 
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Figure 3 | C EDL obtained from direct AC measurement (a), Schematic side view of a device for capacitance measurement and the equivalent circuit. 
Materials and terminal configurations are indicated. R IL and R gl . are the electric resistance in the ionic liquid and graphene, respectively. Cp t is the 
capacitance between Pt wire and the ionic liquid. Actual size of Pt wire (100 urn in diameter) is much larger than the square window (15 um on a side), 
but it is depicted as a thin wire in the figure for clarity. Optical microscope image of a sample is also indicated. Scale bar shows 20 um. (b), V tg dependence 
of C EDL for 6-layer graphene. The gate voltage was swept from negative to positive values. 



This may be one of the reasons why the difference in H\,gllkg takes 
place between two analyses. It is reasonable that Hbgllhg using two- 
terminal mobility enhances from the value (~ 0.5) of Hbglp^g esti- 
mated from Hall bar devices (_R H and a in four-terminal geometry). 
Thus, Hbg/f-kg m two-terminal geometry can be concluded to be close 
to unity as we assumed in the analysis of equation (2). 

Finally, to give physical meaning to the experimental result, we 
calculated the charge distribution in the FLG that forms the EDL 
capacitor, and obtained C EDL as a function of n. We define L7 EDL ( = 
Q*/2Cedl) as the total energy stored in the EDL capacitor where 
charges +Q and — Q are condensed in the ionic liquid and FLG 
boundaries, respectively. As shown in Fig. 4a, the charge — Q is dis- 
tributed over the layers so as to screen penetration of the electric field 
into the interior. The electric field E i between the i-th and i+l-th 
layer becomes smaller than E^ owing to the charge — q i lying on the 
i-th layer, where i = 1, n is numbered from the ionic liquid side. 
Then, C g is obtained by C g = Q 2 /2L/ g , where U g is the static field 
energy due to E^. 



1 



E t {rfdr = 



Q 2 d 



2£il£q 2e, 



where e gr = 5.7 is the relative permittivity of graphite 27 and t = 
0.34 nm is the distance between graphene layers 28 . The first term 
in the right-side expression shows the classical electric field energy 
accumulated in EDL and the second term indicates that in rc-layer 



graphene. Note that a familiar formula C g 
SLG owing to n = 1 



= £/2,£ 0 /d is obtained for 



On the other hand, we can relate C q to the band 
filling energy L/ q 29 , C q = Q 2 /2U q . Originally C q was introduced to 
explain the incomplete shielding of the electric field effect for a two- 
dimensional electron gas 17 . We extend this concept to thicker films 
that can completely screen the electric field by expressing !7 q as a 
summation of the band-filling energy in each layer; 



w=E 



ED(E)dE -- 



E 2e 2 D 

/= 1 



(6) 



where £ Fi is the Fermi energy of the i-th layer. For simplicity, we used 
a constant D(E) to deduce the second equality. As mentioned before, 
C q = e 2 D is obtained for SLG. We minimized Uedl ( = U g + U q ) with 
respect to each qi of the variational parameters and thereby deter- 
mined Ug, L/ q , and L/ EDL . Using these energies, capacitances were 
deduced. C g , C q , and Cedl are plotted as a function of n in Fig. 4b. 
All the capacitances significantly depend on n: C q increases and C g 
decreases with n, and both capacitances saturate to constant values at 
n>4. This saturation behaviour corresponds to the fact that the 
electric field is screened by the charge distributed within 3-4 
layers 2 * 30 in FLG. Note that the charge distribution over 3-4 layers 
in FLG significantly reduces C g because the effective EDL thickness 
expands from d to ~ d + at; a is on the order of unity and determined 
by the effective charge distribution. Therefore, even if an ionic liquid 
with extremely thin d were used, C g could not be increased due to the 
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2 4 6 8 10 
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Figure 4 | Model calculation of C EDL (a) Schematic representation of the charge distribution -qi in FLG, which shields the electric field, 
(b) Calculated capacitances as a function of n. Geometrical capacitance C g , quantum capacitance C q , and EDL capacitance C EDL are indicated by green, 
purple, and black lines, respectively. Details of the calculation can be found in the main text and Supplementary Information. We find that the C EDL is 
limited by C q for n < 4 and by C g for n > 4. 
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additional EDL thickness of at (see Supplementary Fig. S4 online). 
Although layered material is considered to be the most suitable for 
use with an ionic liquid gate, the effect of decreasing C g must be taken 
into account. 

In conclusion, we have clarified the microscopic principles deter- 
mining the capacitance of rc-layer graphene by systematic study. The 
carrier distribution and the band structure in the electrodes of a 
capacitor define C g and C q , respectively, both of which contribute 
to the total capacitance. This study demonstrates that we may not in 
principle achieve the high-density carrier accumulation required for 
E F to reach the van Hove singularity even with ionic liquid gated 
FLG. This result suggests that we need chemical doping 6 ' 31 ' 32 or band 
control 9 to bring the van Hove singularity close to the Fermi energy. 
Our finding for EDL capacitor may establish a nanotechnological 
guideline to improve ultrathin capacitors for charge doping, energy 
storage, and electric power supply. 

Methods 

FLGs were prepared on a Si0 2 /highly doped Si substrate by the micromechanical 
cleavage of Kish graphite (Covalent Materials Co.). The Si0 2 was 300 nm thick and its 
surface was coated with a hydrophobic hexamethyldisilazane (HMDS) layer 33 to 
improve the FET characteristics. The number of graphene layers was distinguished by 
the contrast of the optical microscope image 16,21 " 24 , atomic force microscopy (AFM) and 
Raman spectroscopy 34 " 37 (see Supplementary Fig. SI online). As described in the main 
text, we only used uniform graphene layers without any overlaps of different layers, 
wrinkles and folds in order to obtain a reliable n dependence. FLG devices for two- 
terminal transport measurement were fabricated by photolithography and the vacuum 
evaporation of metal (Cr 5 nm/Au 50 nm) and insulator (LiF 30 nm) to passivate the 
electrodes. FLG devices with a Hall bar structure were prepared using the electron beam 
lithography and oxygen plasma etching. FLG devices for capacitance measurement were 
coated with 800 nm thick photoresist (Tokyo Oka, TSMR8900) instead of LiF, followed 
by photolithography to open a window on the graphene layer. All the measurements 
were carried out in Ar at room temperature. The transport property was measured with 
a semiconductor device analyser (Agilent B1500A). We evaluated the conductivity from 
the slope of J - V curves measured at a drain voltage of ± 1 mV. Except for the transport 
measurement using a Hall bar structure, a gel of bmim[PF 6 ] was used as the top gate 
dielectric (preparation details described elsewhere 38 ). The geometrical capacitance of 
this ionic liquid was estimated to be 9.7 uF/cm 2 by extrapolating the measured capa- 
citance at 20 Hz - 100 kHz to 0 Hz. This value corresponds to d = 0.6 nm, assuming 
En, — 7 !5 . The capacitance measurement was carried out with a precision LCR meter 
(Agilent E4980A) applying an AC voltage of 10 mV at 500 Hz and superimposed DC 
voltage V" tg from — 1 to 1 V. Only the FLG devices with a Hall bar structure were 
measured by applying V*^ with an electrolyte LiClCVPEO. The mixing molar ratio of 
LiC10 4 /PEO is 1/20. The geometrical capacitance of this electrolyte is reported to be 
15 uF/cm 2 39 , which is of the same order as that in bmim[PF 6 ]. 
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